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摘 要 :为 了 考察 高 速 流 场 对 面板 高 频 振 动 的 影响 ,建立 了 高 速 流 场 环境 中 面板 的 能 量 辐射 传递 模 
型 。 针 对 高 速 流 场 环境 中 的 二 维 面板 ,引入 受 高 速 流 场 影响 的 波 数 和 和 群 速度 ,建立 了 能 量 密度 控制 
方程 ,推导 了 高 速 流 场 环 境 中 面板 的 能 量 密度 和 能 量 强度 的 核 函 数 。 实 源 强度 由 导 纳 法 计算 得 到 
的 输入 功率 表示 , 虚 源 强度 通过 边界 的 能 量 平衡 方程 确定 ,根据 串 更 斯 原理 , 板 的 能 量 响 应 由 实 源 
产生 的 直接 场 与 虚 源 产生 的 反射 场 线性 县 加 得 到 。 最 后 ,通过 将 本 研究 所 提 能 量 辐射 传递 法 
(RETM ) 的 计算 结果 与 解析 解 对 比 , 验 证 其 正确 性 ,同时 分 析 了 高 速 流 场 对 板 能 量 响 应 的 影响 。 本 
研究 为 高 速 流 场 环境 下 飞行 器 面板 的 高 频 振动 响应 预测 提供 了 一 种 有 效 的 分 析 方 法 。 
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A study on the energy radiation transfer model for high frequency 
vibration panels in high-speed flow field 
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Abstract :Im order to investigate the influence of high-speed flow field on high-frequency vibration of pan- 
els,a new energy radiation transfer model is established. For the two-dimensional panel in the high-speed 
flow field environment ,the wave number and group velocity affected by the high-speed flow field are intro- 
duced ,the energy density control equation is established , and the kernel functions of the energy density 
and energy intensity of the panel in the high-speed flow field environment are deduced. The real source 
strength is expressed by the input power calculated by the admittance method , and the virtual source 
strengths are determined by the energy balance equation at the boundary. According to Huygens principle ， 
the energy response of the plate is obtained by the linear superposition of the direct field generated by the 
real source and the reflected field generated by the virtual sources. Finally ,by comparing the results calcu- 


lated by the proposed energy radiation transfer method (RETM ) with the analytical solutions ,the correct- 
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ness of the proposed model is verified ,and the influence of high-speed flow field on the energy response of 


panels is analyzed. This research presents an effective analysis method for the high frequency vibration es- 


timation of the flight vehicle panels in high-speed flow field. 


Key words:energy radiative transfer model; high-speed flow field; high frequency vibration; energy re- 


sponse 


随 着 高 速 飞行 器 的 发 展 ,高 速 流 场 环 境 中 面板 
结构 的 高 频 振动 特性 研究 广 受 关注 02 。 高 速 流 场 
环境 中 的 面板 结构 受到 气动 力 载荷 的 激励 ”9 ,引起 
气动 弹性 效应 5 ,从 而 影响 面板 构件 的 振动 特性 。 
因此 ,对 于 高 速 飞 行 中 面板 结构 的 高 频 振动 特性 分 
析 , 必 须要 考虑 气动 压力 载荷 的 影响 。 

实际 上 ,气动 压力 和 其 他 动 载荷 作用 下 面板 的 
振动 响应 分 析 一 直 以 来 广 受 关注 。IBRAHIM 等 
种 用 有 限 元 法 ,采用 一 阶 活塞 气动 理论 分 析 气 动 压 
为 (研究 了 功能 梯度 板 在 流 场 环境 中 的 响应 。 
WANG 等 外 采用 有 限 元 法 和 虚 功 原理 ,分 析 了 超 音 
速 气流 对 金属 板 模 态 相互 作用 和 振动 特性 的 影响 ， 
庆 局 WANG 等 上 又 考虑 了 热 载荷 的 影响 。HAN 
等 中 考虑 板 长 度 方向 均匀 流体 运动 ,推导 了 各 个 方 
国 静波 数 .附加 质量 以 及 等 效 群 速度 ,建立 了 均匀 流 
场 作 用 下 薄板 高 频 振动 响应 的 能 量 流 模型 ,计算 了 
下 极 的 能 量 密度 分 布 ,研究 了 水 流 对 受 集 中 力 激励 
板 的 高 频 振动 特性 的 影响 。WANG 等 "利用 大 拨 
度 板 理论 和 活塞 空气 动力 学 理论 ,建立 了 高 速 流 场 
环境 中 板 的 气动 弹性 运动 方程 ,对 黏 弹性 加 热 板 在 
气 恩 载荷 作用 下 的 非 线 性 颤 振 进 行 了 分 析 。 王 晨 
等 区 | 基于 统计 能 量 分 析 方法 ,分 析 了 不 同 温度 条 件 
下 @@ 行 器 铝 合 金 壁 板 的 高 频 声 振 响 应 。ZHANG 
等 上 9 考虑 了 流体 载荷 引起 的 附加 质量 效应 和 辐射 
阻尼 效应 ,推导 了 流体 载荷 作用 下 平板 的 能 量 控制 
方程 。CHEN 等 5 提出 了 一 种 新 的 能 量 流 模 型 ,用 
于 预测 超 音 速 气流 中 壁 板 的 高 频 振动 响应 ,得 出 结 
果 与 相应 的 精确 解析 解 进行 比较 。 

传统 的 有 限 元 法 和 边界 元 法 被 广泛 使 用 于 低频 
噪声 和 振动 的 预测 中 ,车 将 它们 应 用 于 结构 高 频 振 
动 会 遇 到 两 个 问题 ;由 于 振动 波长 较 小 导致 计 算 
成 本 升 高 ;@ 结 构 不 确定 性 非常 敏感 导致 计算 鲁 棒 
性 不 强 ' 申 。 为 了 克服 传统 确定 性 方法 的 频率 限制 ， 
工程 中 已 发 展 了 多 种 高 频 分 析 方 法 。 其 中 最 成 熟 的 
分 析 方 法 是 统计 能 量 分 析 法 (statistical energy analy- 
sis,SEA)'"1 ,但 SEA 只 能 预测 各 子 系统 平均 能 量 
响应 ,无 法 预测 子 系统 内 部 的 能 量 响 应 分 布 "1。 


而 LE BOT 提出 的 能 量 辐 射 传递 法 '””1, 是 由 热 辐 
射 传递 比拟 而 来 ,属于 几何 声学 范畴 ,可 以 得 到 子 系 
统 内 部 任意 位 置 的 能 量 密度 和 能 量 强度 。 王 幸 '”” 
采用 能 量 辐射 传递 法 分 析 了 梁 和 板 的 高 频 振动 响应 
其 结果 与 波 传播 分 析 法 . 模 态 着 加 法 及 能 量 有 限 元 
法 进行 对 比 ,显示 结果 更 加 准确 。 钟 强 ” 将 RETM 
进行 了 扩展 ,研究 了 各 向 异性 二 维 板 高 频 振动 特性 ， 
并 将 结果 与 解析 解 进行 对 比 。 因 此 ,能 量 辐射 传递 
法 是 一 种 高 效 ,准确 的 高 频 振 动 响应 分 析 方 法 。 

鉴于 现 有 的 能 量 辐射 传递 模型 没有 考虑 气动 弹 
性 效应 ,本 研究 在 LE BOTI*” 1 和 钟 强 ” 的 前 期 工 
作 基 础 上 ,将 高 速 流 场 环境 引入 面板 模型 ,建立 了 高 
速 流 场 环境 中 面板 的 能 量 辐 射 传递 模型 ,同时 分 析 
了 高 速 流 场 对 面板 高 频 振动 响应 特性 的 影响 。 


1 高 速 流 场 环境 中 面板 的 振动 分 析 
模型 


图 1 高 速 流 场 环 境 中 受 简 谐 载荷 激励 的 板 
Fig.1 Plate excited by harmonic load in high-speed flow field 


图 1 为 高 速 流 场 环境 中 受 简 谐 载荷 激励 的 面 
板 , 两 端 简 支 , xo 位 置 受到 横向 简 谐 点 激励 ,产生 弯 
曲 振动 ,其 运动 方程 可 以 表示 为 


4 -2 
De Pl = Fe”6(x — xo) 


(1) 
式 中 ;w 为 板 的 横向 位 移 ; D. = D(1 +jm) 为 复 弯 曲 
刚度 ,这 里 D = E17 为 板 的 弯曲 刚度 , 为 弹性 模 
量 , 1 为 板 横 截 面 的 惯性 矩 , nm 为 阻尼 损耗 系数 
(n 过 1) ; p, 为 板 的 面 密 度 ; Foe”6(x - x) 为 横向 
简 谐 载荷 ,这 里 为 振幅 , w 为 圆 频率 , 6(x - x) 
是 Delta 函数 。 


对 于 高 速 流 场 Mae (V2,5) ,面板 上 下 压力 差 P 
- P, 通过 一 阶 活塞 理论 表示 为 
2 
P-P,= pa ( 
Ma -1 


OWw + Me = 这 aw] 
ox Ma -1 V oa 
(2) 
式 中 : p, 为 空气 密度 ;了 为 气流 速度 ; Ma 为 马赫 数 ， 
且 Ma = Ja , a6 为 音速 。 
将 式 (2) 代 入 式 (1) ,可 以 得 到 高 速 流 场 中 板 的 
弯曲 振动 方程 ” 
Ow Ow 


Ow Ow j 
2 +p, 有 二 大， + C, Fe”6(x — xo) 


D 


(3) 
式 中 : K，= p,V/VMa -1 表示 空气 动力 刚度 ; 
C==sp.V(Ma? - 2)/ (Ma -1)” 表示 空气 动力 
门 波 的 传播 


中 由 于 受到 横向 简 谐 载荷 的 激励 ,在 稳定 状态 下 ， 
合 兆 以 波 的 形式 传播 ,将 w = Ae 代入 式 (3) 
可 得 高 速 流 场 环境 中 板 横向 振动 的 频 散 关系 中, 即 
CN Dk -pw +jKk+ijCw =0 (4) 
求 萎 方程 可 得 4 个 根 ,对 应 一 对 传播 波 和 一 对 修 逝 
波 5 则 式 (4) 的 根 为 "“ 


x ， - ee 
GG Yo ? 


1 


Db. 
oo 4 2 到 
- 天 _ Ps 十 JK,k, ]Csw 
本 四 a D. 是 
4 7 - 
天 Ps + Kk jC,w 
el D. ? 
4 7 
s@ + Kk, jC,w 
天 ， 至 2 (5 ) 


式 中 : 局 表示 右 行 传播 波 波 数 ; h, 表示 左 行 传播 波 
波 数 ; 4. 表示 右 行 俱 逝 波 波 数 ; fo 表示 左 行 做 逝 波 
波 数 。 因 此 ,高 速 流 场 环境 中 板 的 弯曲 振动 位 移 的 
通 解 可 以 表示 为 
w(x,t) = (Ae I + Bei> + Ae 十 ee)eiw 
(6) 

式 中 4 、B 、4. 和 B. 为 待定 系数 ,代入 边界 条 件 可 
以 确定 。 

对 于 右 行 传播 波 = 局 ,+ jia ,代入 式 (5) ,可 
以 表示 为 “ 
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op KR/w), Kk.+Cw 
| De 

(7) 
定义 p! = p. - Kivo 为 右 行 传播 波 的 等 效 密度 ， 
定义 m1 = (一 K+ Csw) /piw 为 气动 压力 产生 的 
附加 阻尼 。 定 义 m，= 7 + 1m 为 右 行 传播 波 的 等 效 
阻尼 。 考 虑 到 mi。 和 为 小 阻尼 , 即 远 小 于 1 的 假 
设 , 所 可 以 表示 为 


1 ~ | (8) 


根据 式 (8) ,可 以 得 到 右 行 传播 波 访 的 实 部 六 和 虚 
部 司 的 表达 式 "” , 即 


4 7 
[9 r 
hi, ~ / 守 ， ki; 一 一 Tk, (9) 


同 理 ,可 以 得 到 左 行 传播 波 h 的 实 部 h, 和 虚 部 hb; 
的 表达 式 '“ , 即 


4 


k,, A 


Pa” 
万 ， 
这 里 , 左 行 传播 波 等 效 密度 p,， = p. + 及 ji[w” ,气动 
压力 附加 阻尼 m,。= (Kk, + Csw)/psw” ,等 效 阻尼 
M2 = 1 + 2a o 
将 上 =- 记 和 有 = 万 分 别 代入 频 散 方程 式 (4) ， 
取 方 程 实 部 ,可 以 得 到 
Di -piw = 0， 
Di -pw =0 (11) 
从 频 散 方程 式 (6) 中 可 以 得 到 , 群 速度 和 相 速 


Vor 
hk; 一 一 hz (10) 


Cgl oh, 本 pi@ ? 
2Dok; 
(12) 
Ok, p20 k, 
2 Dh, 
> hk 四 pr 
Dj 
(CU 2r 
Gs (13) 
hk,, Pp2 


可 以 发 现 ,高 速 流 场 作用 会 影响 面板 波 的 传 


高 速 流 场 作用 中 面板 的 能 量 密度 WW=W, + W,， 
其 中 叹 表示 动能 密度 ,，W, 表示 势能 密度 。 则 时 间 
平均 变量 (于 》、《W,)》 的 表达 式 为 ” 
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Cv) = 了 (2 人 (2 (14) 


式 中 : 〈(*) 表示 时 间 平 均 量 ; " 表示 为 复 共 罗 变 量 。 
因此 ,时 间 平 均 能 量 密度 可 以 表示 为 “ 


(CW) =《W> + (Wi) 一 Pio 


14 |2e26z 


(W) = (Wa) + (Ws) ~3ps0 ?|B |?e2ar 


(15) 
式 中 :下 标 1 号 表示 右 行 传播 波 ; 下 标 2 号 表示 左 行 
传播 波 。 
高 速 流 场 作用 下 板 的 能 量 强度 由 剪 力 产生 的 能 
量 强 度 和 弯 矩 产生 的 能 量 强度 组 成 。 时 间 平均 变量 
( 语 、\ gi J 2 


的 


A 


| 


( 
(加 ) ] ‘16) 
左 行 传播 流 (1,) 可 表 


(Jy) = 


es 
RN 14] 
《| 

DD. (1,) = (I) + (lw) ~ Dok, |Ble™* (17) 
> 根据 式 (15) 和 式 (17), 可 得 能 量 速度 表达 
太仓 

© C1) 
ke (W) 


es (18) 


Dol 


因此 ,能 量 速度 和 群 速度 等 价 , 进 而 可 得 出 能 量 
传递 方程 , 即 
gi =calW), 
(L,) =cp W,) (19) 
1.3 输入 功率 


时 间 平 均 输入 功率 可 以 通过 导 纳 法 


定 59 , 即 


1 1 
P, = 7 Re(F oo) 一 


进行 确 


IF | Re(Y) (20) 


式 中 : vw 表示 激励 点 处 的 速度 ; 了 = vo/P 表示 激励 
点 处 的 导 纳 。 

对 于 高 频 振动 分 析 , 可 以 使 用 无 限 结构 推导 的 
导 纳 近似 有 限 结构 的 导 纳 。 对 于 轴 力 作用 下 无 限 结 


构 二 维 面板 ,激励 点 处 导 纳 了 表示 为 
_@ (ka hk,) (hk, -hk ) — (kk, -ka) (ky -hk,) 
Dk —h) ks -hk) Ck hk) Ck kh,) (hk, -hs) 
(21) 
式 中 ,hk = 一 ki ,ke =k ,Ek =- jk ,ks = -jk o 
将 式 (21) 代 入 式 (20) 得 到 的 Pi, 不 是 有 限 结构 
面板 的 精确 输入 功率 ,而 能 够 表征 精确 输入 功率 频 
域内 的 整体 变化 趋势 。 


2 高 速 流 场 作用 下 二 维 面 板 的 RET™ 
模型 


本 研究 能 量 辐射 传递 模型 研究 的 主要 任务 是 高 
速 流 场 环 境 中 二 维 面板 的 能 量 密度 控制 方程 ,推导 
其 能 量 密度 和 能 量 强 度 核 函 数 ,然后 根据 患 更 斯 原 
理 ,用 实 源 产生 的 直接 场 与 虚 源 产生 的 反射 场 线性 
谷 加 表示 板 的 能 量 响应 。 


六 党 膏 


él 
天 Bd 
a > 
(a) 右 行 传播 波 
A 
OL 
L L+ Be 
a 
(b) 左 行 传播 波 
图 2 高 速 流 场 环 境 中 面板 能 量 平衡 示意 图 


Fig.2 Energy balance diagram of panel in high-speed flow field 
如 图 2 所 示 , 针 对 左右 行 传播 波 微 元 中 的 能 量 
守恒 ,建立 能 量 平衡 关系 式 , 即 


ol 
+ Pi =P.,6(%—%0), 


a 
a 十 Pas 
OX 


式 中 :了 是 x 方向 能 量 强 度 ; Pj, 是 结构 自身 阻尼 产 
生 的 耗 散 功率 ; P;, 是 激励 点 处 的 输入 功率 。 
对 于 振动 频率 为 w ,阻尼 损耗 系数 为 7 的 简 谐 


=P,6(%— xXo) (22) 


弹性 振动 系统 ,其 单位 体积 的 时 间 平 均 耗 散 功 率 与 
时 间 平 均 能 量 密度 成 正比 , 即 
Pisa) = MK Wy, 
(Pi ) = W,) (23) 
将 式 (23 ) 与 式 (19) 代 入 式 (22) ,可 以 得 到 高 速 
流 场 环境 中 二 维 面板 的 能 量 密度 控制 方程 , 即 


ol W» 
ca + mo W) =P.,6(% — xo), 
9 W, » 
co tno W) =P,6(x -x) (24) 


考虑 到 高 频 振 动 响应 分 析 , w 较 大 ,高 速 气流 对 


群 速度 影响 较 小 , 群 速度 可 以 近似 为 
cisco 一 0 =2(0Do]p (25) 
一 在 非 激 励 点 位 置 可 以 表示 为 
2 a 
Ox 
2 +m(W,) =0 (26) 


式 中 mi = w/cs ,ma = mw/c。。 求 解 式 (26) 可 
以 得 到 能 量 密度 表达 式 , 即 

CD (W)> = Ce™! lwo 

(W,) = Ce "1 (27) 


QI 相应 的 能 量 强度 表达 式 为 

~> (7) = Cicse ™! [eso we 

SC (〈 厂 》 = Gere 21| (28 ) 
(6 在 激励 点 x = x 处 ,能 量 密度 是 连续 的 ,能 量 强 


度 泵 连续 ,满足 条 件 


(Wi» | = 《W,) a 
《7 | ,t+ 《7 | 。 = (29) 


=%0 


求 得 CG = Pi,/2c,,C, 汪 Pi/2c, ,进而 能 量 密度 为 


ey 
L 人 in 2c, » 
ee | x-xo | 
CWay = Pe (30) 
从 而 可 得 能 量 密度 和 能 量 强度 的 核 函 数 为 
e™! | x-xo | 
MX > Xo 
Cs 
G(xo,%) = (31) 
ez | x-xo | 
2c. ,XTX < Xo 
El | x-xo | 
2 ew ,> Xo 
H(xo ;%) = | (32) 
-2 x—x0 | 
2 ev X < Xo 
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式 中 e, ,为 xo 到 接收 点 x 的 单位 向 量 。 对 于 无 限 结 
构 二 维 面板 ,激励 点 右 侧 表现 为 实 源 以 右 行 传播 波 
传播 , 左 侧 表现 为 左 行 传播 波 传播 。 

Si 


x=0 X=L 
图 3 高 速 流 场 环 境 中 面板 RETM 模型 示意 图 
Fig.3 RETM model of panel in high-speed flow field 
对 于 有 限 结构 二 维 面板 (图 3) ,考虑 长 度 为 上 的 板 
RETM 模型 示意 图 ,在 x =xo 处 存在 振幅 为 oo 的 实 源 ， 
在 x=0 和 «= 工分 别 存在 振幅 为 os 和 os 的 虚 源 。 
这 里 输入 功率 P。= Si 。 两 侧 虚 源 强度 分 别 用 


os 和 os 表示 。 因 此 ,在 x 点 的 能 量 密度 和 能 量 强度 
可 以 通过 实 源 和 虚 源 的 能 量 和 加 给 出 ,表示 为 
W(x) = aoC(xo,Y) + oaG(A,x) + opG(B,x%), 
I(x) = ooH(xo,x) + oAH(A,x) + osH(B,x) 


(33) 
式 中 oo = P;, 。 能 量 密度 核 函 数 G 和 能 量 强 度 核 函 
数 豆 按照 源 与 测 点 的 位 置 进行 匹配 。 因 此 , 虚 源 or 
表现 为 右 行 传播 波 , 虚 源 os 表现 为 左 行 传播 波 , 面 
板 上 任何 一 点 都 有 右 行 传播 波 和 左 行 传播 波 。 


7 


图 4 高 速 流 场 环 境 中 面板 RETM 模型 中 能 量 的 传播 
Fig.4 Energy propagation in RETM model of panel 


in high-speed flow field 

对 于 能 量 辐射 传递 法 ,要 用 实 源 产生 的 直接 场 
与 虚 源 产生 的 反射 场 线性 厂 加 表示 振动 响应 。 
此 ,需要 通过 边界 处 能 量 平衡 求解 出 边界 处 的 虚 源 
强度 , 即 入 射 能 量 等 于 反射 能 量 ,如 图 4 所 示 。 当 只 
考虑 传播 波 时 RETM 模型 中 能 量 的 传播 ,可 得 
lim(Ts mm) = [h(x = 0) +Ils(x = 0)1m,, 
lm mp) = [L(x = + = 0)1:n (34) 
式 中 :ni 表示 从 4 点 到 B 点 的 单位 向 量 ; z 表示 从 
B 点 到 4 点 的 单位 向 量 。 

将 核 函数 表达 式 (32) 代 入 式 (34) 可 得 


A LP. H(xo,A) + opH(B,A)l1: n,, 


5 = [PH(x,B) +o,H(A,B)1:n, (35) 


02 
SA 这 
联 立 式 (35 ) 可 以 解 得 
已 -2x0 平 emaC-20) —m2l 
OA 一 P;, ] - Pe ? 
—m2x0—m1l -m1(l-x0) 
已 @ 
Oh P;, 1 a @ ml (36) 


将 两 侧 虚 源 强 度 表达 式 (36) 代 入 能 量 密度 和 
能 量 强 度 表达 式 (33) ,可 以 求 得 高 速 流 场 中 二 维 面 
板 上 任意 一 点 的 能 量 密度 和 能 量 强度 。 


3 数值 算 例 分 析 


为 了 验证 本 研究 建立 的 高 速 流 场 中 二 维 面板 的 
能 量 辐 射 模 型 的 正确 性 ,分 析 高 速 流 场 对 二 维 面板 
的 高 频 振动 响应 的 影响 。 研 究 对 象 采用 图 1 所 示 高 
速 滨 场 环境 中 两 端 简 支 的 面板 ,用 能 量 辆 射 法 进行 
高 频 振 动 响应 分 析 , 并 将 得 到 的 能 量 密度 分 布 与 波 
传 轿 法 解析 解 进行 对 比 。 
面板 参数 如 下 :p =2 700 kgjm ,已 =70 GPa, = 
1G® 1 =1 mm,v = 0.3,7 = 0.02, 其 中 p, = 
0C194 75 kg/m’ ,oo = 295. 07 m/s。 高 速 流 场 Ma = 
3 5 身 向 简 谐 力 以 F, =1N 的 振幅 施加 在 板 的 中 点 。 
3 和 高 速 流 场 中 二 维 面板 能 量 密度 
CN 
让 根据 传统 板结 构 富 振 分 析 " ,计算 板 的 临界 马 
赫 数 为 3.2, 所 以 本 研究 中 所 分 析 的 马赫 数 需 要 满 
中 Do <3.2。 


= 393.61 


"之 393.60 


回 393.59 


Sod DD 
Re] Ne] 
id 
tn tn 
的 co 


和 群 速度 /ms 


TO D4 50 
习 5 等 效 群 速度 随 马 赫 数 变化 (f=4 kHz) 
Fig.5 Variation of equivalent group velocity with 
Mach number(f =4 kHz) 

在 本 研究 中 建立 高 速 流 场 环境 中 二 维 面板 的 能 
量 辐 射 传递 模型 时 ,忽视 了 高 速 流 场 对 面板 等 效 群 
速度 的 影响 。 如 图 5 所 示 , 激 励 频 率 4 kHz 时 不 同 
马赫 数 下 面板 的 等 效 群 速度 ,可 以 看 出 左右 行 传播 
波 的 群 速度 受到 高 速 流 场 影响 而 不 同 。 右 行 传播 波 
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的 群 速度 大 于 不 受 流 场 影响 时 面板 的 群 速度 ,并 且 
随 着 马赫 数 的 增 大 而 增 大 ; 左 行 传播 波 的 群 速度 小 
于 不 受 高 速 流 场 影响 时 面板 的 群 速度 ,并 且 随 着 马 
赫 数 的 增 大 而 减 小 。 但 是 ,与 不 受 流 场 影响 时 比较 ， 
高 速 流 场 对 等 效 群 速度 的 影响 较 小 ,不 到 0.01% 。 
图 6 给 出 Ma =3 时 面板 群 速度 随 频 率 变化 图 ,可 以 
看 出 ,高 速 流 场 对 等 效 群 速 度 在 所 有 频段 内 都 影响 
较 小 。 这 说 明 可 以 忽略 高 速 流 场 对 面板 等 效 群 速度 
的 影响 ,近似 为 cu = ce = cs, ,因此 本 研究 推导 的 高 
速 流 场 环境 中 面板 的 能 量 辐射 传递 模型 是 合理 的 。 
1 000 n 

900 : 时 

800 归 

700 
600 

500 
400 
300 

200 

100 


群 速度 /(m *s) 


频率 /kHz 
图 6 等 效 生 速度 随 频率 变化 (Ma =3) 
Fig.6 Equivalence group velocity varies with frequency (Ma =3) 
如 图 7 所 示 ,高 速 流 场 Ma =3 时 面板 的 能 量 密 
度 分 布 图 。 此 处 参照 文献 [14] ,所 示 解 析 解 在 1/3 
倍 频 程 内 取 平 均 , 同 时 也 在 一 个 波长 上 进行 局 部 空 


间 平 均 。 可 以 看 出 ,能 量 辐射 传递 模型 分 析 所 得 的 
能 量 密度 ,除了 激励 点 外 ,其 他 位 置 都 吻合 较 好 。 这 
是 因为 能 量 辐射 传递 模型 在 建立 时 忽视 了 修 逝 波 和 波 
与 波 之 间 干 涉 的 影响 ,而 修 逝 波幅 值 衰减 速度 快 ,在 激 
励 点 附近 一 个 波长 内 几乎 衰减 到 0。 同 时 能 量 辐射 伟 
递 法 分 析 所 得 的 能 量 密度 与 传播 波 非 干 涉 的 能 量 密 度 
吻合 恨 好。 说明 能 量 辐射 传递 模型 所 得 的 能 量 密度 的 
物理 本 质 就 是 传播 波 非 干涉 项 的 时 间 平 均 能 量 密度 。 
6.0 


解析 解 
行 波 非 干涉 解 
RETM 


能 量 密度 /(107J .ma 


0.0 0.2 0.4 0.6 0.8 1.0 
位 置 x// 
图 7 高 速 流 场 环境 中 能 量 密度 分 布 (Ma =3/=4kHz) 
Fig.7 Energy density distribution in high-speed 
flow field environment (Ma=3,f=4 kHz) 


3.2 高速 流 场 对 面板 高 频 振动 的 影响 


图 8 给 出 激励 频率 f=3.0 kHz 时 无 限 结构 面板 
在 马赫 数 等 于 3 的 高 速 流 场 环 境 中 与 无 流 场 环境 中 
的 能 量 密 度 分 布 图 。 对 比 发 现 ,无 流 场 环境 中 无 限 
结构 面板 能 量 密度 分 布 关于 激励 点 位 置 对 称 ,而 在 
高 速 流 场 环境 中 ,面板 的 能 量 密度 关于 激励 点 位 置 
不 对 称 ,并 且 激励 点 左 侧 能 量 密度 衰减 幅度 大 于 右 
侧 衰减 幅度 。 图 9 显示 高 速 流 场 环境 中 左右 行 传播 
波 的 等 效 阻尼 随 马 赫 数 的 变化 图 ,可 以 看 出 右 行 传 
播 波 等 效 阻尼 小 于 左 行 传 播 波 等 效 阻尼 ,再 考虑 到 
无 限 结构 面板 只 有 实 源 产 生 的 直接 场 ,激励 点 右 侧 
只 有 右 行 传播 波 , 左 侧 只 有 左 行 传播 波 。 因 此 ,高速 
流 场 环境 中 无 限 结构 面板 的 能 量 密度 右 行 波 衰减 得 


= =0.5 0.0 0.5 Lg 
位 置 w// 


< 图 8 无 限 结构 面板 能 量 密度 分 布 (f=3 kHz) 


图 
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© Fig.8 Energy density distribution of infinite panel 
= (f=3 kHz) 
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16 18 20 22 24 20 28 30 
马赫 数 
图 9 高 速 流 场 环 境 中 等 效 阻尼 随 马赫 数 变化 
(f=3 kHz,Ma=3) 
Fig.9 Variation of equivalent damping with Mach number 
in high-speed flow field environment (f=3 kHz,Ma =3) 

如 图 10 所 示 ,激励 频率 /=3 kHz 时 高 速 流 场 环 
境 中 面板 和 无 流 场 环境 中 面板 的 能 量 密 度 分 布 图 。 
对 比 发 现 ,高 速 流 场 环境 中 有 限 结构 面板 的 能 量 密 
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度 关 于 激励 点 位 置 不 对 称 ,并 且 相 较 于 图 7, 有 限 结 
构 面 板 受 高 速 流 场 影 响 激励 点 位 置 能 量 密度 不 同 于 
无 流 场 环境 时 。 根 据 式 (36) 可 以 看 出 高 速 流 场 对 
虚 源 强度 产生 的 影响 ,并且 左 右 行 传播 波 等 效 阻尼 
不 同 ,因此 高 速 流 场 环境 中 激励 点 处 实 源 场 与 虚 源 
场 琶 加 后 的 能 量 密度 与 无 流 场 环境 中 不 同 。 


图 10 有 限 结构 卫 
(f=3 kHz) 

Fig. 10 ”Energy density distribution of finite panel 
(f=3 kHz) 


密度 /(10*J，m”) 


= 
转 里 


位 置 x// 


图 11 有 限 结构 面板 能 量 密度 分 布 
(f=15 kHz) 


Fig. 11 Energy density distribution of finite panel 
(f=15 kHz) 

图 11 为 激励 频率 f=15 kHz 时 有 限 结构 面板 在 
马赫 数 等 于 3 的 高 速 流 场 环境 中 和 无 流 场 环境 中 的 
能 量 密度 分 布 图 。 可 以 看 出 ,高速 流 场 环境 中 面板 
能 量 密度 分 布 会 随 着 频率 升 高 而 恢复 对 称 性 。 图 
12 显示 马赫 数 等 于 3 的 高 速 流 场 环 境 中 左右 行 传 
播 波 的 等 效 阻尼 随 频 率 的 变化 图 ,可 以 看 出 ,高 速 流 
场 环 境 中 的 等 效 阻尼 随 频率 升 高 而 更 加 接近 无 流 场 
时 的 结构 阻尼 。 因 此 ,高 速 流 场 会 对 面板 的 能 量 密 
度 分 布 产生 影响 ,而 这 种 影响 会 随 着 频率 的 升 高 而 
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图 12 高 速 流 场 环境 中 等 效 阻尼 随 频率 变化 (Ma =3) 


Fig.12 Variation of equivalent damping with frequency 


in high-speed flow field environment (Ma =3) 


4 三 结 论 


(O 〇 建立 了 高 速 流 场 环境 中 板 高 频 振动 的 能 量 辐 射 
伟 遂 异型。 能 量 密度 控制 方程 由 能 量 平衡 方程 导 
由 引 特别 是 推导 了 与 流 场 相关 的 能 量 密度 和 能 量 强 
度 核 函数 。 推 导 了 高 速 流 场 环 境 中 板 的 近似 导 纳 ， 
i 蜡 输 入 功率 ,用 于 表示 实 源 强度 ,然后 通过 边界 能 
晤 泥 衔 确定 虚 源 强度 ,最 后 采用 能 量 传递 法 用 实 源 
FB 的 直接 场 与 虚 源 产生 的 反射 场 线性 释 加 表示 高 
下 漳 场 环 境 中 板 的 振动 响应 。 所 得 结论 如 下 。 
全 1 ) 通 过 与 精确 解 的 比较 ,发 现 使 用 无 限 结构 导 
缩 主 算 的 输入 功率 的 RETM 结果 与 解析 解 和 行 波 非 
干 器 解 吻合 较 好 ,证 明 本 研究 提出 的 高 速 流 场 环 境 
中 恤 量 辐射 传递 模型 的 有 效 性 。 
-2) 高 速 流 场 环境 中 会 对 左右 行 传播 波 的 等 效 蛆 
尼 产 生 影响 ,进而 影响 左右 行 传播 波 的 能 量 密度 分 
布 。 右 行 传播 波 的 等 效 阻尼 小 于 左 行 传播 波 的 等 效 
阻尼 ,从 而 导致 左 行 传播 波 的 能 量 密度 衰减 更 快 。 
3) 当 流 场 马赫 数 不 变 时 ,等 效 阻尼 随 着 激励 频 
率 增 大 接近 无 流 场 阻 尼 , 即 高 速 流 场 对 面板 振动 响 
应 的 影响 会 随 着 激励 频率 增 大 而 减 小 。 
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